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Abstract 

We consider what constraints unitarity and CPT invariance yield on the strong and electromagnetic phases 
entering K — > tttt decay. In particular, we show that the relative size of the electromagnetically-induced 
changes in the / = and 1 — 2 phase shifts in the two-pion final state do not depend on the explicit coupling 
to the 7r^7r~7 channel. This demonstrates that Watson's theorem can be extended to include the presence of 
electromagnetism. We point out the consequences for the general structure of the K —f -Kir decay amplitudes 
in the presence of isospin violation. 
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1. A detailed understanding of the rich phenomenology of if — > tttt decays has remained elusive despite 
decades of effort. Although progress has been made, the dynamical origin of the A/ — 1/2 rule, as well as 
the strength of CP-violating parameter Re (e'/e), is not yet clear. Another, presumably related, puzzle stems 
from the apparent violation of Watson's final-state theorem. Watson's theorem emerges from unitarity and 
CPT-invariance, in concert with isospin symmetry, and implies that the strong phase in K tttt decay ought 
be given by that of tttt scattering. However, the S-wave tttt phase shift difference So — 62 extracted from the 
K TTTT decay modes, using physical masses in the phase-space integrals, is about 57° whereas its value 
from TTTT scattering data, with the help of chiral perturbation theory and dispersion relations, is about 45° with 
an uncertainty of roughly 10% jl], |j. The assumed equality of these quantities is a consequence of isospin 
symmetry, so that the resolution of the discrepancy has been sought in the computation of isospin- violating 
effects. Isospin violation can be generated by both strong (up-down quark mass difference) and electromagnetic 
(virtual photon) interactions, and its effects have been recently studied in great detail, see, e.g., Refs.|p[|-[p^. 
While many interesting results have been obtained and many others are forthcoming, the gap between the phase- 
shift difference obtained from K tttt decay and tttt scattering has thus far eluded a detailed explanation. In 
the framework of chiral perturbation theory, which is the appropriate theoretical tool in this context, many new 
low-energy constants appear in the most general Lagrangian of Goldstone bosons coupled to virtual photons 
and external sources, making certain numerical predictions difficult. It is thus important to explore whether 
Watson's theorem can be extended in the presence of isospin violation. In Ref. ^ Watson's theorem was shown 
to persist through leading order in the up-down quark mass difference, so that the phase shifts from K tttt 
decay and tttt scattering ought be equal to 0((md — m,i)^). Were the phase shifts from if —> tttt decay and tttt 
scattering equal, the empirical phase-shift discrepancy could nevertheless be resolved, for it could be interpreted 
in terms of an additional amplitude in if ^ tttt decay, of |A/| = 5/2 in character In the isospin-perfect 
limit, the K tttt transition in 0{Gf) can be of |A/| — 1/2 or |A/| = 3/2. In the presence of isospin violation, 
a |A/| — 5/2 transition can be realized from ^ ruu effects in concert with a |AJ| = 3/2 weak transition 
or from electromagnetic effects in concert with a |AJ| = 1/2 weak transition. The empirical enhancement of 
the |A/| = 1/2 weak transition suggests that the latter mechanism is of greater importance. The empirical 
|A7| = 5/2 amplitude required to resolve the phase-shift discrepancy is compatible with that expected from 
electromagnetic effects |8], yet it is significantly larger in magnitude than and of opposite sign to that indicated 
by explicit estimates jH 5|. Moreover, including the estimated phase-shift difference from electromagnetism |5j 
exacerbates this discrepancy. These difficulties prompt the consideration of Watson's theorem in the presence of 
electromagnetism, in order to realize what constraints may exist on the strong and electromagnetically induced 
phase shifts in K ^ tttt decay. This is the aim of the present investigation. It has been triggered by the work of 
Bernstein pT[ |, who considered isospin violation in near-threshold neutral pion photoproduction from protons, 
extending the final-state theorem to the situation of three open channels, in that case 7p, ir^p, and Tr~^n, where 
p{n) denotes the proton (neutron). In a similar fashion, we consider three open channels for the decays, 
which are the two-pion final states with total isospin zero and two, denoted as (7r7r)o and (7r7r)2, respectively, 
as well as the inelastic tt^tt^j channel, whose inclusion is required to render the electromagnetic corrections 
to the K (tttt)/ amplitudes infrared (IR) finite. We will construct a general 4x4 S-matrix appropriate to 
this scenario and derive a set of unitarity constraints from it. Our purpose is not a detailed numerical analysis 
of the various isospin-violating effects, but rather the construction of a theoretical framework which would be 
helpful in constraining such calculations. Nevertheless, we will be able to derive consequences from the unitarity 
constraints which thus far have only appeared indirectly in numerical analyses. 

2. First, we must collect some definitions for the discussion of the K ^ tttt amplitudes — we follow the 
notation and conventions of Ref. |^ and refer the reader to that paper for further details. In the isospin limit 
(m,i = md,e = 0), the decay of a neutral kaon into two pions with isospin / equal to zero or two can be 



parametrized via"^^ 



((7r7r)/|7^w|^°) = AIeM^6I), 

{{TnT)i\Hw\K^) = A*exp(tSi), (1) 

where Tiw is the effective weak Hamiltonian for kaon decay. The amplitude Aj is such that Aj — \Ai\ exp(ji^/), 
with ^7 the weak phase associated with the decay to the final two-pion state of isospin /, and Sj is the phase shift 
corresponding to S-wave tttt scattering of isospin /. In the isospin-symnietric limit, Bose symmetry requires 

'^'^Here A/ is i times Aj defined in the Ref. 
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the pion pair to have / = or / = 2. In that hmit, the S-matrix for strong scattering in the (tttt)/ final state 
is described by a pure phase, 

Here, we have tacitly assumed that at ^/s — M^a the inelasticities from the opening of the four-pion threshold, 
27r — > 47r, are negligible. This is not only a well-known empirical fact but it can also be understood in 

the framework of chiral perturbation theory, for it first occurs in three-loop order (see Ref. |Q, e.g.). 

We now turn to the inclusion of isospin-breaking effects. For the moment we neglect the presence of the 
channel K 7r+7r^7, and we consider merely how isospin violation impacts the tttt subspace. As previously, 
we introduce the channels (7r7r)o and (7r7r)2, where these states are related to the physical basis via^ 



Itt+tt ) cx |(7r7r)o) + -^|(7r7r)2) 

cx l(w)o)- V2|(w)2) . (3) 

Enforcing unitarity and time-reversal invariance, the general S-matrix appropriate to scattering in the two- 
pion subspace — with zero net charge — contains exactly three parameters. Two parameters characterize tttt 
scattering in the isospin-perfect limit, so that the additional parameter must be at least of 0{md — my) or 
0{e), as isospin-breaking is generated by both strong and electromagnetic effects. Our neglect of the 7r"'"7r~7 
channel would be appropriate were we to consider strong-interaction isospin-violating effects only. Let us do 
this and examine the extensions necessary for the treatment of electromagnetism later. Working in analogy to 
the "bar phase shifts" for J ~ S ~ 1 nucleon-nucleon (NN) scattering in the presence of a tensor force , we 
parametrize the S-matrix as 

/ e*^o \ / cos 2k i sin 2k\ / e**'" q \ 
~ V 6*^=^ / V i sin 2k cos 2k / V e'^^ / ' 

where k is the third parameter which is sensitive to isospin-violating effects. In the absence of isospin violation, 
K = 0, and we have Si = Si. G-parity arguments show that strong isospin-violating effects in TTTr-scattering 
are of 0((md - m„)2) 

We parametrize the K tttt amplitudes in the presence of isospin violation via 

{{7nT)i\S\K°) = tAiexpitSi), 

{{7nT)i\S\K^) = tA*j expiiSi) , (5) 

noting that the Si are the strong phases of the K —>■ tttt amplitude and recalling that S = I + iT. Unitarity 
constrains the explicit relation between the Si and the Si, note Ref. If the channel-coupling parameter were 
zero, then Si = Si = Si, and the strong phase in the K —f tttt decay would be that of tttt scattering in the 
isospin-perfect limit. For later use, we introduce the abbreviation 

Ai = Si -Si, 1^0,2, (6) 

so that A/ = for K = 0. 

To end this discussion, we briefly return to the isospin-perfect limit. The 3x3 S-matrix describing the coupling 
of the Kq to the (tttt)/ channels takes the form, where Kq, (7r7r)o, and {1^)2 refer to rows/columns 1,2 and 3, 
in order: 



S = 



( 1 iA^e'^" iAie-'^-" 
iAoe*^" e^''^° I , (7) 



\ iA2 e'"^ 



*^Note that the properly symmetrized state \n 7r+)sym = (Itt^ttj ) + Itt-^ 7r^))/v'2 = v^Itt+tt ). 
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so that ~ {{'^^)o\S\K^) . The ampUtudc Aj contains a non-trivial weak phase. By CPT invariance, 
((7r7r)/|Tt|ifO) = (((7r7r)/|r|i^o))*, so that 5i2 = (is:°|5|(7r7r)o) = ((7r7r)o|S'|ifO). We work in ©(Gf, 6°), though 
we neglect terms of 0{Gf) in our parametrization of the (tttt)! <-> (tttt)// S-matrix. We do this as our interest 
is in the constraints which exist on the T-conscrving phases associated with an amplitude Ai, so that only 
the unitarity constraints emerging from (5^5)12 — (5^ 5) 13 = are of interest. An amplitude Aj is itself of 
0{Gf), so that 0{Gf) contributions to (tttt)/ ^ (tttt)/' scattering play no role in the order of Gf to which 
we work. Thus it is appropriate to neglect 0{Gf) effects in our description of (tttt)/ ^ (tttt)/' scattering; we 
can parametrize this 2x2 matrix by a form which is both unitary and T-conserving. The unitarity constraints 
(iS'''5)i2 = (5''' 5)13 = lead to Si = Sj {I = 0,2), so that the strong phases appearing in K ^ tttt decay are 
exactly those of elastic tttt scattering. This is Watson's theorem in the isospin-perfect world. Equipped with 
these results, we are now in position to generalize this framework to include the 7r+7r^7 final state as well. 

3. We wish to extend Watson's final-state theorem to include both electromagnetic and strong-interaction 
isospin-violating effects in if — > tttt decays. For that, we extend the 3x3 matrix of Eq. to an appropriate 
matrix of larger dimension. Before presenting this extension, let us collect and discuss the assumptions of our 
analysis. 




Figure 1: A schematic illustration of how electromagnetism can generate channel-coupling effects. We let a 
round circle represent a particular final-state channel. The (tttt)/ ^ 7r+7r'7 channel coupling starts in 0{e), 
whereas the (7r7r)o (7r7r)2 channel coupling, mediated by the 7r+7r~7 channel, starts in O(e^). A (7r7r)o <-> (7777)2 
channel coupling can also be mediated by photon exchange, so that the presence of an intermediate 7r+7r~7 
state is not necessary. 



1) We assume that the 7r+7r~7 final state is the only inelastic channel which couples to the two-pion channel 
of given isospin, thus generating transitions of the type (7r7r)o <-> (7777)2 as illustrated in Fig. 1. As we 
have noted, the inelasticities generated by the opening of the four-pion threshold can be safely neglected. 
Considering the empirical branching ratios of Ks to electromagnetic final states, we note that the branch- 



ing ratio for Ks — *■ 77+77^7, presuming photon momenta in excess of 50 MeV/c^, is roughly 2 • 10"'^, 
whereas the next largest measured branching ratio, Ks 77, is roughly a factor of 1000 smaller []l7| . 
Ks — *■ 377 decay is also possible, but is unimportant in this context, as a 377 —> 277 transition with J = 
violates parity and cannot contribute in 0{Gf) to the unitarity relations of interest. Consequently, the 
77+77^7 final state is the only inelastic channel of interest, and the appropriate extension of Eq. (|^) is a 
4x4 matrix. 

2) We work in leading order in the Fermi constant, 0{Gf), and focus on the unitarity constraints in which 
the kaon decay amplitudes, Aj for K — > (7777)/ decay and A^ for K 77"*" 77^7 decay, appear. Thus we 
may choose the 3x3 submatrix without kaons, that is, that of the coupled (7777)0, (7777)2, and 77+77^7 
system, which we term the pion-photon system, to be both unitary and T-reversal-invariant. The 3x3 
matrix contains eighteen parameters; nine parameters are constrained by unitarity, and three more are 
constrained by T invariance, so that the resulting matrix has six non-trivial parameters. Scattering in 
these channels is driven by strong and electromagnetic effects. 

3) Recalling the form of the bar phase shifts in NN scattering we write the 3x3 S-matrix of the pion- 
photon system in an analogous form. A - B ■ A, where ,4 is a diagonal matrix parametrized in terms of the 
phase shifts of 7777 and 77+77^7 scattering, namely Sj and (5^, so that 

/ e'^-' 

A= e''° I , (8) 

I e'^^ 



*®The cut on the photon momentum is required; the tt+tt 7 final state generated by bremsstrahlung from a charged pion is 
infrared divergent. 
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cf. Eq. (Q), and S is a unitary, T- invariant 3x3 matrix containing three parameters. Our form of B 
is inspired by the form of the Kobayashi-Maskawa parametrization of the Cabibbo-Kobayashi-Maskawa 
matrix fl^ . The latter, however, contains 4 parameters, but one can readily define one of the parameters 
in terms of the others to yield a T-invariant matrix. More precisely, we introduce two angles 61,2 and 
one phase S. This assignment will be discussed below. The matrix B is chosen to be: 



B 



/ Cl isi C2 isi S2 

isiC2 cicl + s^e'^ C2S2{ci-e'^) | , (9) 
^isiS2 C2S2{ci-e^^) cisl+c^e^^ 



where we adopt the conventional abbreviations Si — sinOi and Ci = cos0,;. 

4) The form of B in Eq. is compatible with the hierachy of channel couplings. The transitions (7r+7r^7) 
(tttt)/ start at 0{e), whereas the couplings between the two-pion final states (7r7r)o ^ (7r7r)2 are of 0{e^). 
It is seemly, then, that Qi appears as si in the (7r^7r~7) <-> (tttt)/ elements and as ci in the (7r7r)o <-> (7171)2 
elements. Furthermore, at 0((m„ — m^)^), there are no transitions between (tttt)/ and 7r+7r~7, whereas 
there are transitions between (tttt)/ and (7r7r)//. Thus the introduction of the phase S is convenient as this 
can describe the presence of nid ^ w„ effects as distinct from the electromagnetic isospin-violating effects 
characterized by Oi. Of course 5 can contain electromagnetic contributions as well. Note that taking 
01 = 0, (5 = sets all the channel couplings to zero. Moreover, including the parameter 5 as a phase 
means that in the limit in which only 0{S, ©i) terms are kept, all the off-diagonal terms are imaginary. The 
remaining parameter, O2, can be thought of as characterizing the difference of the inelasticity parameters 
in the / = and 1 = 2 channels due to the presence of the third channel. 

5) We work with IR- finite amplitudes throughout. For a detailed discussion of the extraction of the IR- finite 
parts from the full amplitudes, we refer the reader to Ref. 1^. Here, it suffices to say that the potentially 
troublesome contributions can be factored and are thus not relevant to our discussion. 

We can now give the generalized 4x4 matrix, where the if", 7r+7r~7, (7r7r)o, and (7r7r)2 channels are associated 
with rows/columns 1,2,3, and 4, respectively. This assignment is prompted by the foregoing remarks. We thus 
have 



S 



( 1 iA;e*^-' iAle^~^° iA^e''^^ \ 

(10) 



iA^e'^-' ae'^'^y isi C2 e'^^°+^-''> isi S2 e'^^-'+^^'> 
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iAoe'^^" isi C2 e'(^o+^^) {a 4 + sje'^) e^'^" C2 S2{ci - e'^) e 
^^^2 6*^2 isiS2e'(^-'+^^) C2S2(ci - e**)e*(^«+^2) (ci + cle*"^') e^*^^ j 

Based on this matrix, we can now derive the consequences of the extension of Watson's theorem including 
electromagnetism and strong-interaction isospin violation. To the best of our knowledge, this form of the 
unitarity constraints has not appeared previously in the literature. 

Before working out the unitarity constraints realized from this 4x4 matrix, we briefly discuss the relation of 
its 2x2 submatrix for 7r7r scattering to existing parametrizations in the literature. Were the inelastic channel 
not present, three parameters would suffice in characterizing it. We have five parameters, whereas the following 
parametrization for the tttt transition matrix in the isospin basis is proposed in Ref. |^ : 

Only T invariance is imposed, and the "bar" denotes the IR finite amplitudes. The parameter a controls the 
isospin mixing, and the opening of possible new channels, yielding a violation of unitarity in the (7r7r)/ sector, is 
parametrized in terms of the inelasticity parameters 770 and 772 . The assumption of having only one additional 



open channel (7r+7r~7) leads to a correlation between 770 and 772, as seen in Eq. (10). We should also note 
that in the limit ©i = 0, our 2x2 submatrix for the 7r7r system is characterized by four parameters, so that 
one of the parameters is redundant, as the resulting submatrix is both unitary and T-reversal-invariant. The 
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resulting relations between the parameters of this 2x2 matrix and those of Eq. (Q) can be determined, but are 
not transparent. 

4. We proceed to derive the explicit form of the unitarity constraints from Eq. (^0|). Specifically, (5^5)2i = 
(5+5)31 = (5+5)41 = yields 

A-^ = ci e^'^^^-^-') - Ao isi C2 e'(^o~^«+^"'-^"') - A2 is^ s2e*(^2-52+5^-5^) ^ (12) 
= -A^isx C2 e'(*^o-^o+*.-*.) + (ci + si e~'*') e^^^^""^"") + A^ C2 S2 (ci - e''^) e^(^~o-5o+52-5,) ^ 

(13) 

^2 = -A^isi S2 e'Cs^-'^^+~s.-s^^) + C2 S2 (ci - e"'^) e*(^«-*'"+^^-^"^) + ^2 (ci sl + cl e''^) e^'^'^'-^'^ . 

(14) 

For 81 = (5 = we recover ^0 = ^0, ^2 = ^2, and = as we would expect. Recalling Eq. (^ and defining 
we proceed by eliminating A-y exp(iA-y) from Eqs. ( p"3[p^ to recover 



Ao S2 {2i sin Ao + {e-'^ - 1) e*'^") - A2 C2 {2i sin A2 + {e-'^ ~ 1) e'^^) ^ . (15) 
Remarkably these formulae do not depend on 81, and we can factor e^*"^/^ to obtain 

4a.tan8,r^^^^i^i^^ 



Aq Vsin(A2 - (5/2)/ 

We find, as discussed previously, that A2 ^ Aq since the 1 — 2 phase shift is enhanced by a factor of 
A0/A2 ~ 22 1^. Moreover, no terms in 81 appear, so that the explicit coupling to the 7r+7r~7 channel is 
irrelevant to A/. Interestingly, 81 does not enter either Eq. ( p^ or Eq. (p^). This was also observed in Ref. 
but arose from the results of a numerical analysis, whereas here it emerges as a consequence of unitarity. Let 
us make one more comment about Eq. ([l^ ) before proceeding. The right-hand side is explicitly real, whereas 
the left-hand side is not, as the amplitudes Aj carry weak phase information. Requiring the imaginary part 
of the left-hand side of Eq. ( |l6|) to be zero yields the constraint ImA2/Re A2 — ImAo/ReAo = 0, apparently 
suggesting that the CP- violating parameter Re (e'/e) is zero as a consequence of unitarity. However, this is 
not the case: the amplitudes in Aj contain isospin violation as well, and Eq. ([l^ ) becomes indefinite in the 
isospin-perfect limit. If we proceed to examine the relationship between Ai and the amplitudes appearing in 
the general paramctrization of Ref. of the isospin-breaking effects in the K — > tttt amplitudes, we find that 
Re A/ and luiAi arc distinguished by an additional, isospin-violating, CP-conserving function. Physically this 
implies that our paramctrization of the 5^1 and 5ii matrix elements is not sufRciently general, that in the 
presence of isospin violation, there is an additional, path-dependent CP-conserving piece. For our purposes we 
can neglect this additional contribution, though it could well impact the Standard Model prediction of Re (e'/e), 
and we thus proceed to neglect weak phases throughout. Thus we interpret A^ and Aj as KeA^ and as Re A/, 
respectively. 

We can also eliminate A2 from Eqs. ( p2| , p^ to yield 

A^ C2 [cos(A^ + (5/2) - Cl cos(A^ - (5/2)] = Aq si sin(Ao - S/2) , (17) 

and, similarly, Aq from Eqs. ( p^ , p^ , 

A-f S2 [cos(A^ + (5/2) - Cl cos(A^ - S/2)] ^ A2 si sin(A2 - (5/2) . (18) 

It is worth noting that in Eqs. (^^^8|) the parameter 81 does explicitly appear, controlling the relation between 
Aj and Aj. Combining these latter two equations yields Eq. (|l6|), as it should. Since the Aj and A.y are 
complex, our remarks concerning the most general paramctrization of Ai apply to the A^ amplitude as well. 

We also point out that A^ itself is only non-zero in C(e^), S-y being given by the / = l,i = 1 tttt phase shift, 
so that additional expressions are possible. We refrain from reporting these. Nevertheless, we have extended 
Watson's theorem to include the presence of electromagnetism, for the parameters of (tttt)/ and 7r"'"7r~7 scattering 
suffice to relate the electromagnetically generated phases in K ^ tttt. 
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5. In this letter, we have considered the unitarity constraints on the strong and electromagneticaUy induced 
phases in K ~* nn decays. Assuming that the (7r"''7r~7) final state is the only inelastic channel, and working 
in 0{Gf), we have derived the constraints between the decay amplitudes Aj for K — > (tttt)/ decay (for isospin 
/ = 0, 2) and for K tt^tt^j. The corresponding S-matrix is characterized by three mixing parameters; we 
choose two angles &1.2 and one phase 6. The angle Qi is chiefly responsible for the channel couplings (7r+7r^7) ^ 
(tttt)/, whereas S describes the strong isospin violation effects due to the light quark mass difference, as well as 
any "direct" electromagnetic coupling between the (7r7r)o (7r7r)2 states. From the general 4x4 S-matrix for K'^ 
decays into these three final states, cf. Eq.(|o|), one can derive a set of unitarity constraints. Most remarkably, 
it can be shown that the explicit coupling of the (7r"'"7r~7)-channel is irrelevant to the difference between the 
T-conserving tttt phases measured in tttt scattering and extracted from K — > tttt (and K — > 7r+7r~7^^ decays. 
We have also argued that in the presence of isospin violation, the parametrization for the complex-valued decay 
amplitudes Aj ought be modified; we illuminate the sources of isospin breaking which give rise to the general 
parametrization of Rcf. As a next step, it will be interesting to work out the numerical consequences of 
these constraints on the determination of the |A/| —5/2 amplitude in K tttt decay, as well as the impact on 
the Standard Model prediction for the CP-violating parameter Re(e'/e). 
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